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It has been well established that spatially extended, bistable systems that are driven by an 
oscillating field exhibit a nonequilibrium dynamic phase transition (DPT). The DPT occurs when 
the field frequency is on the order of the inverse of an intrinsic lifetime associated with the transitions 
between the two stable states in a static field of the same magnitude as the amplitude of the 
oscillating field. The DPT is continuous and belongs to the same universality class as the equilibrium 
phase transition of the Ising model in zero field [G. Korniss et al., Phys. Rev. E 63, 016120 (2001); 
H. Fujisaka et al., Phys. Rev. E 63, 036109 (2001)]. However, it has previously been claimed that 
the DPT becomes discontinuous at temperatures below a tricritical point [M. Acharyya, Phys. Rev. 
E 59, 218 (1999)]. This claim was based on observations in dynamic Monte Carlo simulations of a 
multipeaked probability density for the dynamic order parameter and negative values of the fourth- 
order cumulant ratio. Both phenomena can be characteristic of discontinuous phase transitions. 
Here we use classical nucleation theory for the decay of metastable phases, together with data from 
large-scale dynamic Monte Carlo simulations of a two-dimensional kinetic Ising ferromagnet, to show 
that these observations in this case are merely finite-size effects. For sufficiently small systems and 
low temperatures, the continuous DPT is replaced, not by a discontinuous phase transition, but 
by a crossover to stochastic resonance. In the infinite-system limit the stochastic-resonance regime 
vanishes, and the continuous DPT should persist for all nonzero temperatures. 

PACS numbers; 64.60.Ht, 75.10.Hk, 64.60.Qb, 05.40.-a 



I. INTRODUCTION 

Metastability and hysteresis are exhibited by numerous natural and artificial systems that are driven away from 
thermodynamic equilibrium by an external "field." In the absence of such a field, and below some critical temperature 
or analogous control parameter, a large class of such systems possess two equivalent, symmetry-broken ordered phases. 
The external field selects one of these ordered phases as the global, stable minimum in the multidimensional free- 
energy landscape. The other ordered phase becomes metastable, separated from the basin of attraction of the stable 
' phase by a free-energy barrier. If the system is initially prepared in the metastable phase and thermal fluctuations 
^ , are present, the system eventually (possibly after an extremely long time) escapes from the metastable free-energy 
I ' well and approaches stable equilibrium. In the present paper we consider the response of a .spatially extended, bistable 
' O [ system driven by an external field which is periodic in time. In particular, we focus on the finite-size effects of the 
^ ■ periodic system response. 

Ferromagnets are perhaps the most commonly known systems that exhibit metastability and hysteresis. In this 
paper we therefore use magnetic language, in which the order parameter is the system magnetization m, and its 
>^ conjugate field is the external magnetic field H. Analogous interpretations, e.g., using the terms "polarization" and 
K> , "electric field" for ferroelectric systems, "coverage" and "electrochemical potential" for adsorbate systems, etc., are 
■ straightforward 

I The dynamic response under an oscillating external field can be viewed as a competition between two time scales: 
the half-period ti/2 of the external field (proportional to the inverse driving frequency) and the average metastable 
lifetime (r) of the system (the mean time spent in the metastable well) after a sudden field reversal. For low driving 
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frequencies, the time-dependent magnetization oscillates about zero in synchrony with the external field (symmetric 
dynamic phase). For high frequencies, however, the magnetization does not have time to switch sign during one 
half-period, and it oscillates about one or the other of its degenerate zero-field values (asymmetric dynamic phase). 
This symmetry breaking and the corresponding dynamic phase transition (DPT) between the symmetric and the 
asymmetric limit cycles of the system magnetization have attracted considerable attention over the last decade. It 
was first observed during numerical integration of the mean-field equation of motion for the magnetization of a 
ferromagnct in an oscillating field [§, Hi- Since then it has been the focus of investigation in numerous Monte Carlo 
simulations of kinetic Ising systems §^1, ||, |, 0, |, g ^ |l| , further mean-field studies §M |, 0, 

and most recently in analytic studies of a bistable time-dependent Ginzburg-Landau (TDGL) model [|17[. The DPT 
may also have been experimentally observed in Co on Cu(OOl) ultrathin magnetic films ^ and recently in 

numerical studies of fully frustrated Josephson-junction arrays and anisotropic Heiscnberg models The results 
of these studies agree that there exists a genuine continuous phase transition between the symmetric and asymmetric 
dynamic phases, at least in some region of the parameter space spanned by temperature, field, and half-period. 
Finite-size scaling studies of data from dynamic Monte Carlo simulations 0, as well as analytic arguments 

have demonstrated that this far-from-equilibrium phase transition belongs to the same universality class as the 
equilibrium Ising model in zero field. This result is consistent with previous symmetry [ p3| and renormalization group 
|P4| arguments ||25|] . 

In their original paper on the DPT in a mean-field model Tome and de Oliveira reported that the continuous 
(second-order) phase transition observed at high temperatures appeared to change at a tricritical point (TCP) to a 
discontinuous (first-order) transition for low temperatures. Such a TCP was also reported in later mean-field work 
1^. However, an analytical and numerical mean-field study by Zimmer |16 makes a strong case that the claims for 
a TCP in the mean-field case is based on a misinterpretation of effects of critical slowing-down at the DPT. Similar 
claims, that in some region of the dynamic phase diagram spatially extended kinetic Ising models exhibit a first-order 
transition and consequently have a TCP separating lines of second-order and first-order dynamic phase transitions, 
have also been made on the basis of dynamic Monte Carlo studies [|[ |2^, |2^. For recent reviews on the DPT, see 
Refs. [HH. 

The purpose of the present paper is to clear up the remaining confusion about the interpretation of simulation 
results for the DPT in spatially extended kinetic Ising models, in particular in the low-temperature regime where a 
first-order transition and a TCP have been claimed to exist 26, Those conclusions were essentially based on 



data for a single system size, and we here demonstrate how proper consideration of rather subtle finite-size effects 
leads to a different picture. The implication of our theoretical arguments and Monte Carlo simulations presented in 
this paper is that in an infinitely large system a continuous DPT should persist down to arbitrarily low temperatures. 
However, in any finite system for sufficiently low frequencies, the DPT gives way to a transient regime of stochastic 
resonance (SR) |28| at a size-dependent temperature. It is this size-dependent crossover temperature, which has 
previously been misinterpreted as a TCP. 

The rest of this paper is organized as follows. In Sec. II we summarize the theoretical framework needed to 
understand the underlying metastable decay mechanisms and their consequences for the DPT. This underscores again 
the importance of the interplay of various time- and length scales in metastable systems ^ |ll|, |3^. In Sec. [II 



we extend our preliminary numerical work [p2[ , supporting our theoretical arguments by large-scale Monte Carlo 
simulations of a two-dimensional kinetic Ising ferromagnet in an oscillating field. Our conclusions are summarized in 
Sec. |v| , and derivations of analytic approximations for quantities of interest in the stochastic-resonance regime are 
given in Appendix ^ and Appendix 



II. METASTABLE DECAY MODES AND PERIODIC RESPONSE IN FINITE SYSTEMS 



The appropriate dynamic order parameter in the presence of an oscillating external field is the period-averaged 
magnetization, Q= 2t\/2 /"^(O*^^ lH- It takes a nonzero value in the asymmetric dynamic phase, while it vanishes in 
the symmetric phase. The transition occurs when the half-period ti/2 and the underlying metastable lifetime (r(T, H)) 
become comparable. The metastable lifetime depends on the temperature T and the field amplitude H. For sufficiently 
large systems (see quantitative statements below) the system escapes from the metastable phase through the nucleation 
of many droplets [multi-droplet (MD) regime p9| , ^]. Consequently, the time-dependent system magnetization is 
self-averaging. If (t(T, H))<^ti/2, the magnetization follows the external field in each half-period. The system relaxes 
to a symmetric limit cycle, and the order-parameter probability density P{Q) is sharply peaked about Q—Q. On the 
other hand, for {t(T, H))^ti/2 the magnetization does not have enough time to switch within a single half-period, 
and it relaxes to an asymmetric limit cycle (with occasional switches between the two equivalent asymmetric dynamic 
phases). Consequently, P{Q) becomes bimodal with sharp peaks near Q=ztl. This breaking of the symmetry of 
the limit cycle and the associated DPT have been carefully analyzed R, O, |4| with finite-size scaling techniques 
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borrowed from equilibrium critical phenomena p3| , In terms of the dimensionless half-period, Q=ti/2/ {t{T^ H)), 
the DPT occurs at a critical value <dc^O{l). The finite-size scaling analysis of the Monte Carlo data also indicates 
that this far-from-equilibrium phase transition belongs to the same universality class as the equilibrium Ising model 
in zero field. Supporting these numerical studies, recent analytic results within a coarse-grained TDGL model 
indicate that the behavior of the stochastic variable Q is governed by the effective Hamiltonian HcS—aQ^ + Q^, where 
aoc{Q — Oc). According to standard arguments from the theory of critical phenomena, this result leads directly to 
the conclusion that the DPT belongs to the universality class of the zero-field Ising model in equilibrium in 
agreement with the simulation results. 

For any finite system, however, the metastable decay mode changes to the nucleation and growth of a single droplet 



at sufficiently low temperatures [single-droplet (SD) regime ||29|, 30 ]. Due to the stochastic nature of the nucleation 
of a single droplet, the corresponding response in the presence of an oscillating field is different: the system exhibits 
lH, |ll|, ^ stochastic resonance ||2^ . 

The crossover from the underlying MD to SD decay can be understood by using standard nucleation theory applied 
to the ferromagnetic kinetic Ising model (in general dimension d and with ferromagnetic coupling constant J) |o| . 
Below the critical temperature, following a single, instantaneous field reversal, the average time between nucleation 
events (the nucleation time) in a system of linear size L is obtained as 

U=[L''liT,H)]"\ (1) 

where I(T, H) is the temperature- and field-dependent nucleation rate per unit volume. It can be expressed in terms 
of the free energy F{T, H) of the critical droplet as 

/(T,iJ)-C(T,iJ)-ie-^(^'^)/^, (2) 

where F[T^ H) and the prefactor C(r, H) can be obtained from nucleation theory to various degrees of approximation, 
depending on T and H ]30| . The temperature T is given in energy units by setting the Boltzmann constant k-Q—l in 
Eq. (^. The other characteristic time scale is the growth time ig. It is defined as the time it takes for a supercritical 
droplet to grow to fill half the system volume. Assuming a time-independent radial growth velocity v{T,H), 

kiL, T, H) = [2f^^(y)]i/d^(j.^^) ' (3) 

where VLd{T) is a dimension- and temperature-dependent shape factor (between tt and 4 for d^2 and between 47r/3 
and 8 for ci=3) 

For tn<^tg many droplets nucleate while those nucleated shortly after the field reversal are still growing. In this, 
the MD, regime the lifetime is independent of the system size. An order-of-magnitude estimate of the metastable 
lifetime {t{T,H)) can be obtained |Q by equating the volume Rq reached by a droplet that grows for this amount 
of time, Rq(x{v{t))'^ , with the volume inside which on average a single nucleation event occurs during the same time, 
Rq(x{{t)I)~^ (Rq can be regarded as the typical droplet separation). The result is (r) oc {v'^I)~^^^'^~^^\ 

For tn^tg the first droplet to nucleate eventually fills the system on its own. In this regime the lifetime depends 
strongly on the system size and approximately equals t^ [Eq. (||)] . 

The above two decay modes characterize the MD and SD regime, respectively. The crossover between the two 
regimes defines the dynamic spinodal (DSP) |^0| and can be estimated by equating and tg. (In terms of the 
underlying length scales, this corresponds to the situation in which the mean droplet separation Rq — v{t) becomes 
comparable to the system size L.) This yields an implicit equation for the temperature corresponding to the DSP as 
a function L and H, 

rp ^ F{Tusp,H) 

"^"^ id+l)lnL-ln{C{Tusp,H)[2n4Tusp)V^''viTusp,H)} - ^' 
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For d—2, analytic approximations |36) are known for v(T,H) [^, Q2(T), and C{T,H 
estimate for Tdsp as a function of L at = 2.0J, obtained by numerically solving EqT^L 
curve in Fig. n^. The resulting curve Tusp{L) turns out to be quite insensitive for the precise values used in the 
approximatiorTpq] . In the large-system limit Tdsp approaches zero logarithmically with increasing L. The existence 
of the DSP implies that in the presence of an oscillating field, reducing T at fixed half-period ti/2 and field amplitude 
H results in drastically different behaviors for "small" and "large" systems. 

Upon reducing the temperature for sufficiently large L, the L- independent metastable lifetime {t{T,H)) becomes 
comparable to the fixed half-period at a temperature above the one at which the underlying decay mode would 
cross over to the SD mode. The horizontal line in Fig. |l| corresponds to the temperature where {t{T,H)) w ti/2 = 
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FIG. 1: Metastable decay modes and crossovers in the dynamic phase diagram for the two-dimensional kinetic Ising ferromagnet 
in a square- wave oscillating field with amplitude H—2.0J. The dashed curve is the dynamic spinodal (DSP) TDSp(i), separating 
the underlying multi-droplet (MD) from the single-droplet (SD) regime. The solid curves correspond to Tx{L) [Eq. the 
crossover to the dynamic "frozen" state for various half-periods ti/2, indicated by arrows [in units of Monte Carlo steps per 
spin (MCSS)]. The horizontal dotted line indicates the temperature T=0.216J, at which the dynamic phase transition (DPT) 
occurs in the large system-size limit for fi/2=500 MCSS. Along this horizontal line, for L>C'(10^) the metastable lifetime is 
system-size independent, and its value is comparable to the fixed half-period. 



500 MCSS. Near this line, the system exhibits a continuous DPT Further lowering the temperature has no 

effect on the behavior of the dynamic order parameter Q, since the system is already locked into one of its symmetry- 
broken dynamic phases, where P{Q) is sharply peaked near ±1. From Fig. ^ it is clear that, in order to observe the 
DPT in a finite system for a given value of ti/2, one has to employ a sufficiently large system [e.g., L>C'(10'^) for 
^1/2=500 MCSS or L^O{10^) for ii/2=50 MCSS]. In general, the larger ti/2 is, the larger systems one needs in order 
to observe the L— »oo behavior. This example illustrates that it is essential to keep in mind the finite-size implications 
of the crossover between the MD and SD decay mechanisms. 

Upon reducing the temperature for smaller L, the crossover from the MD to the SD regime occurs at a temperature 
well above that at which the infinite-system DPT would occur. In the SD regime the switching is stochastic [ pO[ . 
The system will then exhibit stochastic resonance, and the order-parameter distribution P{Q) will indeed possess 
multiple peaks, leading to a negative fourth-order cumulant [^. Lowering the temperature still further reduces the 
probability that the magnetization switches during a finite number of periods, and the system becomes effectively 
"frozen." The crossover to this frozen phase is approximately marked by the curve along which the probability is 1/2 
that the magnetization does not switch during a half-period in which it is opposite to the direction of the field. This 
crossover curve can be estimated by noting that SD switching is brought about by the nucleation of a single droplet, 
which is a Poisson process. Thus, the probability that the magnetization does not change sign during a half-period 
in which it started off opposite to the field is 



p , r exp{-[ii/2-tg(i, T,i7)]A„(L,T,i7)} for tg<ti/2 

Setting Pnot = 1/2 leads again to an implicit equation for the corresponding crossover temperature. 
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(5) 



(6) 



Estimates of at H=2.0J for three different values of ti/2, obtained by numerically solving Eq. (^) |Q, are shown 
in Fig. as solid curves. Since Pnot changes rapidly between zero and unity as T is reduced, the specific cut-off value 
of Pnot — 1/2, used to define Tx here, is not essential. For each value of ti/2, the curves representing Tdsp and Tx 
form the border of a wedge-shaped region in which SR is observed (see Fig. m. This is the regime where the dynamic 
order parameter is indeed characterized by a multi-peaked probability density and a negative fourth-order cumulant 
1^. However, the above theoretical arguments imply that the SR behavior is a finite-size effect. We therefore conclude 
that the otherwise sound simulation results that Acharyya obtained in the stochastic regime , were misinterpreted 
by him as signs of a first-order DPT. 
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To summarize the periodic response for the ferromagnetic Ising model, for an infinite (or sufficiently large) system, 
the system undergoes a continuous DPT when the the half-period becomes comparable with the metastable lifetime. 
Crossing the dynamic phase boundary (Fig. [|) by changing the temperature, for fixed and low frequencies (long half- 
periods), the metastable lifetime becomes comparable to the half-period at appropriately low temperatures, thus the 
DPT occurs at a low temperature (horizontal dotted line in Fig. For finite and too small systems, however, the 
underlying decay mode crosses over to the SD regime before the DPT occurs (crossing the dashed curve from above 
in Fig. 1^). Then the system exhibits SR in the wedge-shape region, until practically no switching occurs during any 
finite observation time (frozen state). As we show next, careful analysis of simulation results for systems of different 
sizes reveals that the signatures attributed to a first-order DPT in Ref. ||^ indeed disappear in the L^oo limit, as 
predicted by the theoretical arguments presented above. 



III. SIMULATION RESULTS AND FINITE-SIZE EFFECTS 



To model spatially extended bistable systems we performed dynamic Monte Carlo simulations of a two-dimensional 
kinetic Ising ferromagnet below its equilibrium critical temperature. This simple model has, for example, been shown 
to be appropriate for describing magnetization dynamics in highly anisotropic single-domain nanoparticles and uniaxial 



thin films ||18|, 19, gO|, |40|. Despite its simplicity, it is believed to capture the generic features of periodically driven, 
spatially extended bistable systems. The system, which is defined on a two-dimensional square lattice of linear size 
L, is described by the Hamiltonian 

n = -jY,s,s,-H{t)Y,s,, (7) 

where Si—±.l is the state of the ith spin, J>0 is the ferromagnetic coupling constant, j) runs over all nearest- 
neighbor pairs, and H[t) is an oscillating, spatially uniform applied field. We use a square- wave field with amplitude H . 
This has obvious computational advantages over a sinusoidal field, while it docs not change the universal characteristics 
of the system response [ p^ . The dynamic used is the single-spin- flip Glauber algorithm with updates at randomly 



m respo 

chosen sites p3| , At temperature T, each attempted spin flip from Si to ~Si is accepted with probability 

-AB./T 

^(-' - = TT-A^ ' (8) 

where /S.Ei is the energy change that would result from the accepted flip. We give the temperature T in energy units 
by setting the Boltzmann constant k-Q=l in Eq . (B ). For the largest system sizes (L>1024) we implemented a scalable 



massively parallel version of this dynamic [^2|, [43[ , first proposed by Lubachevsky 1 44 1 . 
The dynamic order parameter |^ is the period-averaged magnetization 

Q = <f m{t)dt , (9) 

where m(t) — L^^ ^ iSi{t). The beginning of the period is chosen at a time when H{t) changes sign. In particular, we 
compute both types of period averages, staring at the instant when H{t) changes from +H to —H, and also starting 
when it changes from —H to +H. Both observations are included in the order-parameter histograms with equal 
weight. This simplest form of phase averaging is sufficient to produce a symmetric distribution for Q, in particular in 
the stochastic regime. 

Large-scale simulations and finite-size scaling studies of the DPT have been recently performed with both sinusoidal 
10, |l2| and square- wave Q fields. The results imply that the system undergoes a continuous phase transition as the 
half-period becomes comparable to the average metastable lifetime {t{T,H)). Recall that the lifetime becomes 
independent of the system size for large systems. The critical exponents for the dynamic order parameter and its 
fluctuations at the DPT are consistent with those of the two-dimensional equilibrium Ising transition [ij, ^ . 

In those studies, the temperature and the field amplitude were held fixed, resulting in a fixed lifetime (t(T, H)). The 
DPT was approached by tuning the half-period ti^2 of the oscillating field so that it became comparable to {t{T, H)). 
An advantage of this approach is that if the smallest system is already in the MD regime, all the larger ones are, as 
well. Thus, one does not have to deal with subtle crossovers corresponding to the different underlying decay modes 
(SD vs MD). In the present study, we keep ti/2 fixed and tune the metastable lifetime by varying the temperature 
T. The motivation for this is to closely parallel the study by Acharyya and to show that ignoring the finite-size 
effects and the resulting crossovers, one can easily misinterpret the stochastic-resonance behavior in the stochastic 
regime as indicating a first-order transition. 
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FIG. 2: Magnetization time series m{t) (solid curves) and normalized applied field H{t)/H (dashed curves) with H=2.0J and 
ii/2 = 50 MCSS, shown for a "small" system with L=16 at different temperatures, (a) T—0.8J, dynamically disordered phase, 
(b) T=QAJ, stochastic resonance, (c) r=0.35J, dynamically "frozen" state. 



Tracing the magnetization time series m{t) already reflects the major qualitative differences between the responses 
for "small" and "large" systems, as shown in Figs. ^ and respectively. For sufficiently small L, as the temperature 
is reduced, the system enters the stochastic regime, characterized by occasional random switches [Fig. ^(b)], before 
becoming completely "frozen" [Fig. ||(c)]. For large L, the system undergoes a DPT characterized by the slow 
wandering of the period-averaged magnetization [Fig. ^(b)] on its way to performing an asymmetric limit cycle in the 
dynamically ordered phase at still lower temperatures [Fig. ^(c)]. 

We performed simulations on system sizes ranging from L=16 to 2048, choosing various field amplitudes H and 
half-periods ti/2 that were kept fixed while the temperature T was varied. The time unit used is one Monte Carlo 
step per spin [MCSS], i.e., one random "sweep" of the LxL lattice. The lengths of the runs were 10^ full periods of 
the oscillating field for ti/2= 500 MCSS and 10^ full periods for all the other values of <i/2- Measuring the period- 
averaged magnetization Q after each half-period, we constructed averages of the norm of the order parameter, {\Q\)l- 
We further calculated the fourth-order cumulant ratio. 



which typically provides a strong indication of the nature of any underlying phase transition [p3| . For a continuous 
transition, Ul changes monotonically from to 2/3 as one tunes the system from the symmetric (disordered) to the 
symmetry-broken (ordered) phase. On the other hand, for a first-order transition, Ul develops a minimum, whose 
location corresponds to the transition point. We also constructed histograms of Q, representing the order-parameter 
distribution P{Q). In the stochastic regime, we in addition measured the residence times tj- and constructed their 
probability distribution Pr(ti.), the residence-time distribution (rtd). Here t^ is defined as the time spent in one of 
the two "wells" of the underlying system free energy between two consecutive switching events. It is measured as the 
time elapsed between two consecutive zero-crossings of m(t). The residence times and their distribution are useful to 
characterize the system in the stochastic SD regime. 

Keeping H and ti/2 fixed, we performed simulations measuring {\Q\) and Ul as functions of T for a series of system 
sizes. Here we present the results for four different pairs of values of H and ti/2- Figures |^, ||, |^, and ^ show the 
results for H=2.0J and ti/2=50 MCSS, H=2.0J and ti/2=500 MCSS, H=1.8J and ti/2=20 MCSS, and for H=3.0J 
and ti/2=20 MCSS, respectively. As indicated by these figures, the findings in all these cases are qualitatively the 
same. For the purpose of discussion, we use Fig. ^, corresponding to H=2.0J and ii/2=50 MCSS. 



A. Small systems 



Even for relatively small systems, L=16-128, at sufficiently high temperatures the underlying metastable decay 
mode is MD, as illustrated in the phase diagram in Fig. ^(a). Then the system magnetization m{t) follows a symmetric 
limit cycle [see Fig. ||(a)] . Consequently, the order-parameter distribution P{Q) is sharply peaked about zero [Fig ||(a)] , 
yielding ([Q[)l~0 up to finite-size effects, as shown in Fig. ||(b). Correspondingly, the fourth-order cumulant Ul is 
close to 0, as shown in Fig ^(c). 

As the temperature is lowered beneath Tbspl^), the underlying decay mode crosses over to the SD regime, and the 
magnetization switching becomes stochastic as shown in Fig. 0(b). The order-parameter distribution then has three 
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FIG. 3: Magnetization time series m{t) (solid curves) and normalized applied field H{t)/H (dashed curves) with H=2.0J and 
ii/2 = 50 MCSS, shown for a "large" system with L=180 at different temperatures, (a) T—0.5J, dynamically disordered phase, 
(b) r=0.375J near the DPT. (c) T=0.34J dynamically ordered phase. 
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FIG. 4: Dependence of the system response on T and L for H—2.0J and fi/2=50 MCSS. (a) Metastable dynamic phase diagram 
analogous to Fig. ^. The different curves have the same interpretations as in that figure. The horizontal line corresponds to 
r=0.349J, where in the MD regime {t{T,H)) ti/2 = 50 MCSS. (b) The dynamic order parameter (|Q|), shown vs T for L 
between 16 and 1024. (c) The fourth-order cumulant ratio Ul, shown vs T for L between 16 and 1024. Note that the dip to 
negative values disappears as L increases beyond 128. 
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FIG. 5: Dependence of the system response on T and L for H—2.0J and fi/2=500 MCSS. (a) Metastable dynamic phase diagram 
analogous to Fig. ^. The different curves have the same interpretations as in that figure. The horizontal line corresponds to 
r=0.216J, where in the MD regime {t{T,H)) ti/2 = 500 MCSS. (b) The dynamic order parameter (|Q|), shown vs T for L 
between 64 and 2048. (c) The fourth-order cumulant ratio Ul, shown vs T for L between 64 and 2048. Note that the dip to 
negative values gradually disappears with increasing L. 
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FIG. 6: Dependence of the system response on T and L for H—1.8J and fi/2=20 MCSS. (a) Metastable dynamic phase diagram 
analogous to Fig. ^. The different curves have the same interpretations as in that figure. The horizontal line corresponds to 
r=0.592J, where in the MD regime {t{T,H)) ^ ti/2 = 20 MCSS. (b) The dynamic order parameter (|Q|), shown vs T for L 
between 16 and 1024. (c) The fourth-order cumulant ratio Ul, shown vs T for L between 16 and 1024. Note that the dip to 
negative values disappears as L increases beyond 64. 
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FIG. 7: Dependence of the system response on T and L for H—3.0J and fi/2=20 MCSS. (a) Metastable dynamic phase diagram 
analogous to Fig. ^. The different curves have the same interpretations as in that figure. The horizontal line corresponds to 
r=0.193J, where in the MD regime {t{T,H)) ti/2 = 20 MCSS. (b) The dynamic order parameter (|Q|), shown vs T for L 
between 16 and 1024. (c) The fourth-order cumulant ratio Ul, shown vs T for L between 16 and 1024. Note that the dip to 
negative values disappears with increasing L. 



peaks: two extremely narrow peaks near Q—±l and a rather wide one centered at Q—0 [see Fig. ||(b) and (c)]. The 
peaks near ±1 represent the periods during which the magnetization does not switch, while the peak centered on zero 
represents the periods during which it switches at least once. The large width of this central peak is the result of the 
square-wave shape of the applied field, which results in an exponential probability density for the switching process. 
A sinusoidal field would yield a distribution more sharply peaked about zero, since, in that case, the switching almost 
always occurs when the external field assumes its maximum magnitude |ll|, The generic feature in the 

stochastic regime, regardless of the shape of the driving field, is the multiple-peak structure. For the square- wave field 
used in this paper one can obtain (see Appendix |a|) an analytic approximation for P{Q) in the regime where ii/2S>ig 
[Fig. |l^(a)] 

PiQ) = ^S{Q + 1) + |e-«l«l + ^S{Q - 1) . (11) 

Equation ( pT| ) for P(Q) is compared with simulation data in Fig. |o|(b) for a system with L=32 at H=2.0J and 
T=0.34J, for which (t(T, iJ))i=233 MCSS. (The subscript L in {t{T,H))l is included as a reminder that the 
metastable lifetime depends on L in the SD regime.) The half-period is ii/2=500 MCSS ^ tg. This comparison 
contains no fitting parameters: the average metastable lifetime of the underlying metastable decay, {t{T, H))l, was 
measured in single field-reversal simulations, and its value was used to determine the scaled half-period, 6=^1/2/ {t). 
Instead of the delta-functions with amplitude e^®/2, the finite value e~®/ {2AQ) was used to make the correspondence 
with the finite bin-size AQ, employed to build the histogram for Q. In the stochastic regime the cumulant Ul becomes 
negative, reaching a minimum at some temperature, as shown for the smaller values of L in Fig. ^c). 
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FIG. 8: Histograms representing in a small system (I/=16 with H=2.0J and ti/2=50 MCSS) for different temperatures. 
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FIG. 9: Histograms representing P{Q) in a large system (L=180 with H—2.0J and ti/2=50 MCSS) for different temperatures 



A quantity often used to detect stochastic resonance, is the rtd, P,{t,) [0 H, |l|, H, |6), the probabihty density 
for the residence times between zero-crossings of the magnetization. This quantity is shown in Fig. |l^ for a system 
with L=32 at H=2.0J and r=0.34J, for which {t{T, H)}l=233 MCSS. The half-period is ti/2=500 MCSS > tg. 
The data are shown together with an analytic approximation (ii/2^ig), which is derived in Appendix ^ 



-710 



(r) 1 - e 



-e 




^ - 2(n - l)ej if 2{n - l)ti/2 < t, < {2n - l)h/2 
2ne - ^1 if {2n - l)ti/2 < t, < 2nt,/2 



(12) 



The analytic form contains as a parameter the average lifetime {t{T, H))l, which was measured in independent field- 
reversal simulations. Thus, as for the order-parameter distribution above, no fitting parameters are involved in the 
comparison between the simulation data and the analytical form. The generic feature of the rtd is the structure of the 
peaks, which are centered at odd multiples of the half-period (i.e., at (2n — l)ii/2, n = 1,2, . .X, with exponentially 
decaying heights, as has also been observed in simulations with a sinusoidally varying field |ll|, ^l|. This behavior 
is characteristic of systems undergoing stochastic resonance 

mil,!!!- 

In the phase diagram shown in Fig. |^(a), the weak- noise stochastic-resonance |^ behavior described here occurs 
in the wedge-shaped region between the two crossover curves, rDSp(-^) and Tx{L), as already discussed in Sec. II. 
Analogous behavior was discussed in detail in Ref. |^ for the case of a sinusoidal driving field. 

As T is lowered further for the small system, the average metastable lifetime quickly increases, and the probability 
of not switching during a half-period, Pnotiti/2'i L,T, H), approaches unity. As discussed in Sec. II, the crossover 
temperature Tx{L) corresponds to P^ot = 1/2. Consequently, switching events become rare, and the central peak in 
P{Q) essentially disappears, leaving only the two sharp peaks near ±1 [see Fig. ||(d)]. At the same time, Ul again 
becomes positive [see Fig. ^(c)]. Significantly below Tx (i), switching will never be observed during a finite number of 
periods, and the system is completely "frozen" into one of its two metastable wells [see Fig. ^(c)]. This yields {\Q\)^l 
and Ul~2/3, as shown in Fig. ^(b) and (c), respectively. 



B. Large systems 



For larger systems [L>C'(10^) for ti/2—50 MCSS], at high temperatures the system is deeply in the MD regime 
[Fig^(a)], where the lifetime is independent of L. Here, {t{T,H)) is significantly smaller than tif2- The limit cycle of 
the magnetization is symmetric [Fig. ||(a)], and the fluctuations in Q are Gaussian and centered about zero [Fig. ||(a) 
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FIG. 10: Order-parameter statistics for H=2.QJ, T=0.34J, and L=32, which yield (r(r, ))l=233 MCSS, with 
^1/2=500 MCSS. These are the same parameters as for Fig. |l^. For these parameters the system is in the stochastic regime, 
and tg is neghgible compared to and {t{T,H))l- (a) A short segment of the magnetization time series. The interpretation 
of the line types is the same as in Fig. ^. (b) Comparison of the simulated (histogram) and the analytic [solid curve, Eq. (pT[)] 
order-parameter distributions. 
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FIG. 11: Comparison of the simulated and the analytic [Eq. (p^)] residence-time distributions, Pr(tr) in the stochastic regime 
where < 1^/2- Shown for a small system with L=32 at H=2.QJ and r=0.34J, where {r(r, H)) l=233 MCSS. The half-period. 



^1/2=500 MCSS, is much longer than the growth time, (a) Shown on a linear scale, (b) The same as 
scale to emphasize the smaller peaks and the exponential dependence of the peak heights on tr. 



using a linear-log 



and (b)]. In this regime, both (|Q|)l~0 and ?7l~0, up to finite-size effects [Fig. ^b) and (c)]. Upon lowering the 
temperature, the underlying decay mode remains MD, but the lifetime increases and eventually becomes comparable 
with the half-period. This happens well before T reaches TDSp(i)- When {t{T, H)) becomes approximately equal to 
ti/2, the system undergoes a genuine continuous phase transition, the DPT ^]. The system magnetization 

performs a slow "wandering" motion [Fig, ^(b)] , and the distribution for Q widens significantly [Fig. ^(c)]. Below the 
transition, P{Q) becomes bimodal [Fig. |^(d)], {\Q\)L^O{l)jFig. ||(b)], and Ul approaches its ordered-phase value 
2/3 without exhibiting negative values or a minimum [Fig. W(c) . Also, Ul for different large values of L intersect 
at the temperature corresponding to the DPT [Figs, ^c), ^c), |(c), and|^(c)], as expected for a continuous phase 
transition [Q . Detailed finite-size-scaling analysis of Monte Carlo simulations for systems that are large enough that 
the underlying metastable decay mode is MD, are found in Refs. @, |l|,|l4). 

As T is reduced further, the underlying decay mode crosses over to the SD regime at Ty)sp{L). This leads to 
extremely large metastable average lifetimes, such that Pnot approaches unity. However, this has no effect on the ob- 
servables: below the DPT, the system is already performing asymmetric limit cycles, confined to one of its metastable 
wells [Fig.|(c)]. 
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C. Comparison 



From the above discussion for small and large systems, it is clear that the qualitative behavior observed as the 
temperature is varied depends strongly on the field amplitude, the half-period, and the system size. For example, for 
H=2.0J and ii/2=50 MCSS one must employ L>C'(10^) in order for the underlying metastable decay mode to become 
MD, so that the DPT is observed [Fig. gc)]. For the same field amplitude and ii/2=500 MCSS one needs L>O{10^) 
to achieve the same effect [Fig. |^(c)]. For larger and larger half-periods, the "infinite" -system DPT {{T{T,H))K.t 1/2) 
occurs at lower and lower temperatures. However, at these low temperatures, it takes exponentially large systems to 
be in the MD regime, as seen from Eq. (||). By increasing ti/2^ one therefore quickly reaches the limit of any available 
computational resources. 



IV. SUMMARY AND CONCLUSION 



We addressed the finite-size effects of the periodic response of spatially extended bistable systems by studying the 
two-dimensional kinetic Ising model in an oscillating external field. The intimate connection between the underly- 
ing metastable decay modes and the periodic response of spatially extended bistable systems has been stated and 
demonstrated several times [|], |ll|, 12[ 0. On the other hand, it has been claimed repeatedly [^, ^ ^ that 
for large enough periods (low- frequency regime), the dynamic phase transition (DPT) becomes first-order. In the 
present study we focused explicitly on demonstrating that any signatures resembling a first-order transition at lower 
temperatures are merely finite-size effects that disappear as L is increased sufficiently. 

First, we reviewed the basics of the well-known theory of homogeneous nucleation. Understanding the relevant time- 
and length scales and the various decay modes [multi-droplet (MD) and single-droplet (SD)] in metastable decay, one 
can estimate the important system-size dependent crossovers for the periodic response. Next we presented extensive 
new simulation results indicating that no first-order transition exists for any frequency, and consequently, there can be 
no tri-critical point separating lines of first-order and continuous dynamic phase transitions. The behavior, correctly 
observed in Ref. [^ but misinterpreted as indicating the existence of a first-order DPT, is due to the stochastic nature 
of the underlying single-droplet metastable decay. In this regime the system exhibits stochastic resonance. However, 
this behavior does not survive in the large-system limit with fixed field amplitude. 
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APPENDIX A: ANALYTIC APPROXIMATION FOR THE ORDER-PARAMETER DISTRIBUTION IN 

THE STOCHASTIC REGIME 



In the stochastic regime, for large ti/2, we can neglect the growth time after the nucleation of the critical droplet 
and approximate the switching process by instantaneous switching of m{t) between ±1 [Fig. [lo|(a)]. Then, knowing 
that the nucleation of a critical droplet is a Poisson process (i.e., the probability density of the switching time is 
exponential) with rate (t)~^, we can calculate the probability density function (pdf) for the dynamic order parameter 
P{Q)- It is important to recall that in the definition of Q we included both types of averaging: averaging over a 
period when the driving field starts with negative value and also when it starts with a positive value in the first 
half-period. Here we show the calculation of the former case P-{Q). The calculation for the latter is identical and 
at the end one simply has to add them together with weight 1/2, resembling the way the histogram was collected: 
P(Q) = [P_(Q)+P+(Q)]/2.^ 

When the driving field is negative in the first half-period, there are five mutually distinct scenarios as illustrated 
on the schematic plots in Fig. |l^(a-e). In cases (a), (b), and (c) the value of the magnetization m{t) is +1 at the 
beginning of the period, while in cases (d) and (e) it is —1. Since we are interested in the stationary distribution 
of Q, first we have to find the stationary probabilities p'^ (p^) that the magnetization has the value -fl (—1) at 
the beginning of a period. After a quick look at the five cases above [Fig. |lj(a-e)], one can write down a set of 
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FIG. 12: Schematic plots for constructing the probability density P-(Q) when the driving field is negative in the first half- 
period. Dashed and solid lines represent the driving field and the magnetization, respectively. The switching times t\ and 
t2, possibly occurring in the first and second half-period, respectively, are measured from the beginning of their respective 
half-periods. 



discrete-time "time-evolution" equations (from one period to the next) 



p 



n+l 



Pn+i = (l-Oe-^ + e-^p- (Al) 

for p+ and p^, the probabilities that the magnetization is +1 and —1 at the beginning of the nth period, respectively. 
In Eq. (Al) we used the definition <d—tii2/{T). From these equations one can easily obtain the stationary-state 



(fixed-point) values p+ =1/(1 -I- e~^) and p^=e~®/(l -I- e~®). Now for each case (a-e), we consider the conditional 
probability density of Q (conditional on the value of the magnetization at the beginning of the period). In case 
(a) Q—1 with probability e~®, i.e., that the magnetization does not switch in the first half-period, yielding a delta- 
function contribution e^^5{Q — 1) to the full pdf. In case (b) the magnetization switches at ti in the first half- 
period and docs not switch back in the second, resulting in Q={ti — ti/2) /ti/2- The contribution to the pdf is 
{5{Q — {ti — ti/2)/ti/2))ti, where {■ ■ ■)ti is an averaging over the exponentially distributed switching time ti. 
In case (c) the magnetization switches twice, at ti in the first half-period and at ^2 in the second one, resulting in 
Q=[ti — t2)/ii/2- The contribution to the pdf is {5{Q — {ti — <2)/ii/2))ti,t2J where (. . .)tiM is an averaging over the 
exponentially distributed switching times ti and t2, occurring independently in the first and second half-period (and 
measured from the beginning of their respective half-periods). In case (d) the initial value of the magnetization is — 1 
and it does not switch in the second half-period, yielding e^^5{Q + 1). In case (e) the magnetization switches once 
at t2 in the second half-period resulting in Q=— ^2/^1/2 and a contribution {5{Q -\- i2/^i/2))t2 to the pdf. Combining 
the above conditional pdfs with the probabilities of the corresponding initial values of magnetization, one obtains 

P_(g) = p+ {e'^<5(g - 1) + e-«(5(Q - {t, - h,2)lh/2))u + {KQ - {ti - t2)/h/2))t„t,} (A2) 
{e-'-'SiQ + 1) + {S{Q + t2 Ai/2))t2 } • 
Carrying out the averages above using the exponential pdfs for ti and t2 we arrive at 

^-(Q) = {e-®^(Q - 1) + i?(-Q)ee-«(2-|Qi) ^ I (^g-eiQl _ g-e(2-|Q|)j | ^^3) 



{e-«,5(g + i) + if(-Q)ee-«'i«'} , 



where H{x) is the Heaviside step-function. An identical calculation for the pdf of Q for periods starting with positive 
value of the driving field yields 

P+(g) = -£l^{e-^S{Q~l) + H{Q)ee-^\Q\} (A4) 

+ _!_ |e-e<5(Q + 1) + H(Q)ee-«(2-|QI) + | (,-e|Q| _ ,-e(2-|QI)) | . 

The symmetrized ( "phase-averaged" ) dynamic order-parameter becomes considerably simpler and easier to compare 
with measured histograms 

PiQ) = \ [P+(Q) + P_(Q)] = '-^5{Q + 1) + + ^<5(Q ~ 1) . (A5) 
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APPENDIX B: RESIDENCE-TIME DISTRIBUTION AND ITS ANALYTIC APPROXIMATION IN THE 

STOCHASTIC REGIME 



In the stochastic-resonance hmit, where (r) is not much smaller than ti/2, while both are much larger than ig, we 
can obtain an analytic form for the rtd Pr{tr). The derivation follows that given for a sinusoidally oscillating field 
in the Appendix of Ref. However, the present case is simpler since the probability that a switching event has 

not occurred within a time t after the field has changed its sign to become opposite to the magnetization direction, 
Pnot{t; L,T, H), is a simple exponential, exp{—t/{T{T,H))L) [see Eq. As a consequence, all the integrals that 

have to be evaluated numerically in the sinusoidal case, can here be trivially calculated analytically. Provided the 
magnetization switched in a period (say period n=l) at time ti (measured from the instant the driving field changed 
sign) , the next magnetization switching occurring in the nth period at t2 (also measured from the instant the driving 
field changed sign) results in a residence time — {2n — l)ti/2 — ti + t2, where ti and t2 are exponentially distributed 
variables. The formal expression for the rtd then can be written as 



g-("-i)e 

n—1 



(Bl) 



Carrying out the averages above yields, after some rearrangement. 



P,(t,) = 



(r)l 




^ - 2(n - 1)6 if 2(n - l)ii/2 < t, < {2n - l)ti/2 



if (2n - 1)^1/2 <tr< 2nti/2 



(B2) 



where n—1, 2, . 



S. W. Sides, P. A. Rikvold, and M. A. Novotny, Phys. Rev. E 59, 2710 (1999) and references therein. 

T. Tome and M. J. de Oliveira, Phys. Rev. A 41, 4251 (1990). 

J. F. F. Mendes and J. S. Lage, J. Stat. Phys. 64, 653 (1991). 

W. S. Lo and R. A. Pelcovits, Phys. Rev. A 42, 7471 (1990). 

M. Acharyya and B.K. Chakrabarti, Phys. Rev. B 52, 6550 (1995). 

M. Acharyya, Phys. Rev. E 56, 1234 (1997). 

M. Acharyya, Phys. Rev. E 56, 2407 (1997). 

M. Acharyya, Phys. Rev. E 58, 179 (1998). 

M. Acharyya, Phys. Rev. E 59, 218 (1999). 

S. W. Sides, R. A. Ramos, P. A. Rikvold, and M. A. Novotny, J. Appl. Phys. 79, 6482 (1996). 
S. W. Sides, R. A. Ramos, P. A. Rikvold, and M. A. Novotny, J. Appl. Phys. 81, 5597 (1997). 
S. W. Sides, P. A. Rikvold, and M. A. Novotny, Phys. Rev. Lett. 81 834 (1998). 
G. M. Buendia and E. Machado, Phys. Rev. B 61, 14686 (2000). 

G. Korniss, C. J. White, P. A. Rikvold, and M. A. Novotny, Phys. Rev. E 63, 016120 (2001). 

G. M. Buendi'a and E. Machado, Phys. Rev. E 58, 1260 (1998). 
M. F. Zimmer, Phys. Rev. E 47, 3950 (1993). 

H. Fujisaka, H. Tutu, and P. A. Rikvold, Phys. Rev. E 63, 036109 (2001); Erratum, Phys. Rev. E 63 059903 (2001). 
Y.L. He and G.-C. Wang, Phys. Rev. Lett. 70, 2336 (1993). 

Q. Jiang, H.-N. Yang, and G.-C. Wang, Phys. Rev. B 52, 14911, (1995); J. Appl. Phys. 79, 5122 (1996). 

J. S. Suen and J. L. Erskine, Phys. Rev. Lett. 78 3567 (1997); J. S. Suen, M. H. Lee, G. Teeter, and J. L. Erskine, Phys. 

Rev. B 59, 4249 (1999). 

G. S. Jeon, H. J. Kim, M. Y. Choi, B. J. Kim, and P. Minnhagen, Phys. Rev. B 65 184510 (2002). 



H. Jang and M. J. Crimson, Phys. Rev. E 63, 066119 (2001); H. Jang, M. J. Crimson, and C. K. HaU, e-print cond- 



mat/0205070 



G. Grinstein, C. Jayaprakash, and Y. He, Phys. Rev. Lett. 55, 2527 (1985). 
K. E. Bassler and B. Schmittmann, Phys. Rev. Lett. 73, 3343 (1994). 

The direct correspondence between the DPT and other non-equilibrium phase transitions considered by Refs. ^] 
is somewhat less obvious in that no microscopic dynamics is defined here for the dynamic order parameter (the period- 
averaged magnetization). Further, a priori it could not be known whether or not the period-averaging in time introduces 
long-range effective interactions for the local dynamic order parameter. 
[26] M. Acharyya and B. K. Chakrabarti in Annual Reviews of Computational Physics I, edited by D. Stauffer (World Scientific, 

Singapore, 1994), p. 107. 
[27] B. Chakrabarti and M. Acharyya, Rev. Mod. Phys. 71, 847 (1999). 

[28] L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni, Rev. Mod. Phys. 70, 223 (1998). 



13 



[29] H. Tomita and S. Miyashita, Phys. Rev. B 46, 8886 (1992). 

[30] P. A. Rikvold, H. Tomita, S. Miyashita, and S. W. Sides, Pliys. Rev. E 49, 5080 (1994); H. L. Richards, S. W. Sides, M. A. 
Novotny, and P. A. Rikvold, J. Magn. Magn. Mater. 150, 37 (1995); P. A. Rikvold and B. M. Gorman, in Annual Reviews 
of Computational Physics I, edited by D. StaufTer (World Scientific, Singapore, 1994), p. 149. 

[31] S. W. Sides, P. A. Rikvold, and M. A. Novotny, Phys. Rev. E 57, 6512 (1998). 

[32] G. Korniss, P. A. Rikvold, and M. A. Novotny, in Computer Simulation Studies in Condensed Matter Physics XIV, edited 
by D. P. Landau, S. P. Lewis, and H.-B. Schiittler, Springer Proceedings in Physics Vol. 89, (Springer, Berlin, 2002), in 
press. 

[33] K. Binder and D. W. Heermann, Monte Carlo Simulation in Statistical Physics. An Introduction, 3rd edn. (Springer, 
Berlin, 1997). 

[34] K. Binder, in Fmite-Size Scaling and Numerical Simulation of Statistical Systems, edited by V. Privman (World Scientific, 
Singapore, 1990), p. 173. 

[35] N. Goldenfeld, Lectures on the Phase Transition and the Renormalization Group (Perseus Publishing, Reading, MA, 1992). 

[36] In our numerical solution of Eqs. (^) and (^ we employed the low-temperature approximation for the classical theory 
of nucleation and growth where the shape-, prefactor, and the free energy of the critical droplet are independent of the 
temperature |38|. In particular, for the shape factor we used $^2 ~ 4, corresponding to a square droplet. For the the prefactor 
we used C{T,H)^3/8, 1.36, and 5/4 for H£{J,2J), H=2J, and He{2J,4J), respectively For the free energy of the 
critical droplet we used F{T, H)^16J~6H , 4J, and 8J-2H for H€{J,2J], H^2J, and _H'g[2J,4J), respectively [Q . The 
velocity v{T,H) was obtained from Eqs. (13) and (21) of Ref. Q 

[37] P. A. Rikvold and M. Kolesik, J. Stat. Phys. 100 377 (2000). 

[38] M. A. Novotny, in Computer Simulation Studies in Condensed Matter Physics IX, edited by D. P. Landau, K. K. Mon, and 
H.-B. Schiittler (Springer, Berlin, 1997); M. A. Novotny, in Computer Simulation Studies in Co ndensed Matter Phy sics 



XV, edited by D. P. Landau, S. P. Lewis, and H.-B . Schiittler (Springer, Berlin, in press), e-print cond- mat/0204098 
[39] A. Bovier and F. Manzo, e-print |con"d-mat/010737^ (2001). 
[40] P. A. Rikvold, M. A. Novotny, M. Kolesik, and H. L. Richards, in Dynamical Properties of Unconventional Magnetic 

Systems, edited by A. T. Skjeltorp and D. Sherrington (Kluwer, Dordrecht, 1998), p. 307; M. Kolesik, M. A. Novotny, and 

P. A. Rikvold, Phys. Rev. B 56, 11791 (1997). 
[41] R. J. Glauber, J. Math. Phys. 4, 294 (1963). 

[42] G. Korniss, M. A. Novotny, and P. A. Rikvold, J. Comput. Phys. 153, 488 (1999). 

[43] G. Korniss, Z. Toroczkai, M. A. Novotny, and P. A. Rikvold, Phys. Rev. Lett. 84, 1351 (2000). 

[44] B. D. Lubachevsky, Complex Systems 1, 1099 (1987); J. Comput. Phys. 75, 10 (1996). 

[45] P. Jung, Phys. Rep. 234, 175 (1993). 

[46] L. Gammaitoni, F. Marchesoni, and S. Santucci, Phys. Rev. Lett. 74, 1052 (1995). 



14 



